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This paper considers factors that contribute to poor identification of unsteady 
aerodynamics from wind tunnel data for an airliner configuration. One approach to 
modeling a wing-tail configuration is considered and applied to both steady and large- 
amplitude forced pitch oscillation wind tunnel data taken over a wide range of angles of 
attack but a limited range of amplitude and frequencies. The identified models fit the 
measured data well but in some cases with inaccurate parameters. Only limited conclusions 
can be drawn from analysis of the current data set until further experiments can be 
performed to resolve the identification issues. The analysis of measured and simulated data 
provides some insights and guidance on how an effective experiment may be designed for 
wing-tail configurations with nonlinear unsteady aerodynamics. 
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Nomenclature 

non-dimensional pitching moment coefficient 
non-dimensional normal force coefficient 
mean aerodynamic chord, m 
frequency, Hz 

distance between aircraft c.g. and a.c. of the tail, m 

pitch rate, rad/sec 

reference area, m : 

time, sec 

velocity, m/s 

wing center-of-pressure position, fraction of c 

angle of attack, rad 

dynamic angle of attack, rad 

mean angle of attack, rad 

amplitude of angle of attack, rad 

downwash angle, rad 

tail setting angle, rad 

estimated standard error 

time constants for normal force, center of pressure, and downwash angle 
residuals 

angular frequency, rad/sec 
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subscripts 

A = amplitude 

t = tail 


superscripts 

B,T,W = body, tail, wing components 
BT, BW = body-tail, body-wing components 


Aerodynamic Derivatives 

= where a = N or m 
“ da 

c dC 

C„ = — where a = N or m 

? 2V da 

abbreviations 

c.g. = center of gravity 

a.c. = aerodynamic center 


I. Introduction 

O ne of the first attempts to formulate aerodynamic model equations with unsteady terms based on wind tunnel 
data was reported in Ref. 1 and later in Ref. 2. The experimental data were in the form of forced oscillations in 
pitch with fixed amplitude and varied mean angle of attack and frequency. The test article was a simple delta wing. 
Later in the same decade, several reports and papers appeared with experiments including oscillations in roll and 
yaw, and with varied amplitudes 3,4 . Other planforms were also tested such as the F-16XL aircraft with cranked-delta 
wing 5 and the X-3 1 aircraft with delta wing and canard 6 . 

Recently there has been increased interest in identification (model structure determination and parameter 
estimation) of a transport aircraft represented by a classical wing-tail configuration. The theoretical background for 
model formulation was given by Jones and Fehlner 7,8 and Tobak 9 . Later, Klein 10 developed linear aerodynamic 
equations for planar motion of an aircraft with a horizontal tail. Unsteady effects in these equations were expressed 
by linear indicial functions for wing and tail aerodynamics and for the downwash angle. 

A different model of the wing-tail combination was proposed by Khrabrov 11 et. al. They assumed that the aircraft 
model can be represented by body, wing, and tail components and that the aerodynamic forces and moments are in 
general nonlinear and unsteady. The resulting model included algebraic and differential equations with parameters 
dependent upon the angle of attack. The model was used in the analysis of steady and unsteady wind tunnel data. 
The identified model was presented in several graphs comparing measured and estimated time histories. 
Unfortunately, no numerical values of estimated parameters and their accuracies were given. 

Wind tunnel data and postulated model equations of Ref. 1 1 were used again in model identification and the 
results are summarized in Ref. 12. The identified models for different configurations mostly fit the measured data 
quite well. Some of the estimated parameters, however, were estimated with low accuracy resulting in non-physical 
values. There were no a priori values of the parameters and no additional measured data which could be used for 
model validation. 

The purpose of this paper is to re-examine the results of previous data analysis addressing mainly measured data 
inconsistency, low accuracy and non-physical values of parameter estimates, and model adequacy. This will provide 
the groundwork and direction to investigate appropriate experiment designs that can allow identification of wing-tail 
configurations that exhibit nonlinear unsteady dynamics. In addition, a proposal for future experiment design for 
obtaining wind tunnel data with high accuracy and high information content will be outlined. The experiment should 
provide steady data for all configurations, small amplitude oscillatory data for estimating damping parameters, large 
amplitude oscillatory data for model identification, and specific data for checking model adequacy and prediction 
capability. 


II. Measured Data 

The measured data were obtained from static and dynamic testing of an airliner in the Central Aero- 
Hydrodynamic Institute (TsAGI) low subsonic wind tunnel". The steady data included the normal force and the 
pitching moment of four different configurations, i.e., body alone (B), body-tail (BT), body-wing (BW) and the 
complete model (BWT). These measurements cover the angles of attack from -10° up to 30°. 

The dynamic data were obtained from a forced oscillatory motion of the same model in the same tunnel using the 
same test rig. The experiment was executed at different initial values of angle of attack, do, frequencies, and 
amplitudes, a A . For aircraft model identification, large amplitude (a a > 5°) time histories of C v and C„, for three 
configurations BT, BW, and BWT, at frequencies approximately equal to 0.4Hz, 0.8Hz, and 1.3Hz, were used. For 
large oscillation amplitudes of 15°, a 0 was equal to 5° and 15°, and for oscillation amplitudes of 10°, a 0 was equal 
to 5°, 10°, and 20°. Only one cycle of oscillation was given in each time history. 
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III. Models for Steady Aerodynamics 

The aerodynamic forces and moments of the complete aircraft are assumed to be the sum of separate 
contributions from the body, tail, and wing 11 . Then the corresponding models for the steady case have the form 

C N {a) = C B {a) + C^{a) + S t C T N {a t ) (1) 

C m (a) = C B (a) + C% 0 +C%(a)x w (a)-l t S t C T N (a t ) (2) 

where a t =a + tj>-s{a) (3) 


S t = S t / S , i t = 1 1 1 c , and x W is the distance between the wing center of pressure and aircraft c.g. as a fraction of 
mean aerodynamic chord. The static values of C N and C,„ for a complete model and its components are shown in Fig. 
1 and 2. 



alpha, deg 


Figure 1. Normal-force coefficient for four 
configurations (Ref. 11). 



alpha, deg 

Figure 2. Pitching-moment coefficient for four 
configurations (Ref. 11). 


In these plots B, BW, BT, and BWT indicate the values for the body, body-wing, body-tail, and the complete model. 
The values of Cy and Cy in Eq. (1) and Eq. (2) were obtained from the relations 


_ s~ iBW f^B /a\ 

vv 

C T N =C BT -C B (5) 

In Eqs. (1-3), is equal to the wing pitching moment when Cy is equal to zero, </) is the tail plane setting and e is 
the downwash angle at the tail. From the geometry of the model S t = 0.26, l t = 4.62 and (j) = -10° and from C^ n the 
value of CZ Q = -0.13. 

The variable x W (a) can be computed by removing terms associated with the tail in Eq. (2) and using the 
corresponding BW measurements. The downwash angle follows from solving Eqs. (1-3) for a t at selected a. The 

two variables, x H (a) and e(a) , are plotted in Figs. 3 and 4 and compared with results from Ref. 11. The 
differences between the two pairs are small. 
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Figure 3. Distance x w estimated from steady wind tunnel data. 
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Figure 4. Downwash angle at the tail estimated from steady wind tunnel data. 
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Figure 5. Tail length determined from aircraft 
geometry and from tail coefficient measurements. 

In pitching-moment Eq. (2) the last term { -l t S t C^{a t ) } represents the tail contribution, C„(a t ) . In order to 
determine the best representation for this term, two evaluations of the consistency between the measured normal 
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T — — T 

force and pitching moment were made. First, both C m {a t ) and - l t S t C N {a t ) terms were used to compute the tail 
length as 




C 

' — n 


s,c t n 


( 6 ) 


and plotted against a in Fig. 5. Changes in i t with respect to a should be relatively small. Large values of 

— T 

computed i t around a = 15° are mainly caused by C N approaching zero, see Fig 1. For a > 20°, the differences 
between l t computed and its nominal value are larger than expected, indicating one contributor to data 

T T 

inconsistency may result from inaccuracies of measured C m and C N . 

The second check of data consistency is based on a comparison of normal force and pitching moment of the 
complete aircraft with that computed from separate measurements of the body, body-tail, and body-wing. For 
complete consistency, the following relation applies. 


r BWT _ r B, r T r W 
'-'a — + '-'a 


( 7 ) 


where subscript a is either N or m. Using Eqs. (4) and (5) the relation in Eq. (7) can be written as 


r BWT _ r BT r BW _ r B 
'-'a '-'a + '-'a 


( 8 ) 


A check of the pitching moment equation can also be made by replacing Cj n in Eq. (7) by - > g* v * n g 


r BWT _ r BW _ 7 a r T 


( 9 ) 


Figure 6 shows the comparison of measured C^ WT with that computed using Eq. (8) for both normal force and 
pitching moment and Eq. (9) for pitching moment. Normal force comparisons demonstrate very good agreement 
between measured and computed data. Based on comparisons shown in Figs. 5 and 6, one can expect better accuracy 
T T 

of measured C N than C m . Consequently, Eq. (2) is used in this study as stated. 




Figure 6. Measured and computed coefficients 
for complete aircraft . 
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IV. Analysis of Large Amplitude Oscillatory Data 

The aerodynamic model equations for modeling and data analysis are given in Refs. 11 and 12. They are repeated 
here for convenience. 


C N (t;a) = C^(co;a) + (t;a)+ S t C^(t;a t ) + ^-C Nq (co;a)a(t) 


(10) 


C ffl (t;«) = C ; fK«) + C 0 


+ CJT (t; a)x w (t; a) - l t S t C T N ( t ; a t ) + C mq (oo; a )«(0 


(11) 


where 


= a + </> + a d - s(t;a ) 



r ,v( Q ')Qvx = Cj^ s (oo;a) 
z x (a)x w +x w =x fV (oo;a ) 
r £ (a)^ + ^ = ^(oo;«) 

The normal force coefficient due to the wing is partitioned as 

_ r w , r w 

L N “ ^NA + L jVS 


( 12 ) 

(13) 

(14) 

(15) 

(16) 

(17) 


where C^ A is the linear part (attached flow) and C^ s is the nonlinear part (separated flow). Symbol oo indicates a 
variable evaluated in steady flow. Unknown parameters in the model are the time constants, t n , t x , t £ , and two 
damping terms, C N (oo ;a) and C m ^ (oo ;a) . In general, each of these parameters can be a function of a. Whenever 

necessary these equations were modified to allow estimation with better accuracy. 

A. BT Configuration 

The model equations for the BT configuration follow from Eqs (10-13) as 

C N (t) = C B N (a) + S t C T N {a t ) + C T Nq (a)^ (18) 

C m it) = C B m (a) - l t S t C T N {a t ) + C T mq (a)^ (19) 

a, =a + <!) + a d . (20) 

where a d reflects a tail angle-of-attack increment due to rotation given by Eq. (13). 

In these equations the only unknowns are the damping parameters, C T Nq and C T mq . Using the data for a 0 = 5°, a 0 

= 5°, and/= 0.4 Hz, and assuming that the parameters are not a dependent, the following least squares estimates 
and their standard errors were obtained as 

C T Nq = 6.5±1.7andC^ =-21±30 
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with the coefficient of determination, R ? , close to zero. These results suggest inaccurate measured data and/or 
modeling errors in Eqs. (18 or 19). 

To find the reason for such high parameter inaccuracy and low R 2 values, the residuals v N and v m were 

computed from Eqs. (18-19), respectively and examined. Plots of the residuals against time and a are presented in 
Figs. 7 and 8. 




as a function of time. as a function of angle of attack. 

From these plots the dependence of residuals on a is evident. For that reason additional terms were added to the 
model equations yielding 


C,v (0 = C,v ( a ) + S,C t n (a , ) + A Na a + C T Nq (a) ^ 


(21) 


C m {t) = Cl (a) - i t S t C T N {a t ) + A,„ a a + C T mq (a) ^ 
The new parameter estimates in Eqs. (21-22) were obtained as 


C T Nq =6.5 ±0.73 

C T mq=- 21 ±2.0 

A n =0.053 ±0.0022 

iy a 

A =-1.006 ±0.0059 

m a 

R 2 =0.84 

R 2 =0.99 


( 22 ) 


A comparison of results indicates the following: 

a) The mean values of damping parameters remain the same because in harmonic motion a and q are 
orthogonal. 

b) The use of a new model resulted in a dramatic improvement of parameter accuracy and model adequacy. 

c) Estimated bias terms, A and A , are close to a linear approximation of the slope for v N and v m , see 

Fig. 8. 

For the remaining discussion, the estimates will be based on Eqs. (21-22). 


7 

American Institute of Aeronautics and Astronautics 



a, deg 

Figure 9. Measured normal-force and pitching moment 
coefficients for three frequencies at a 0 = 5° and a A = 15°. 

For an assessment of frequency dependence, the measured C N and C„, for three frequencies were plotted in Fig. 
9. In this figure a minor effect of frequency can be seen. Its significance was evaluated by a comparison of two sets 
of estimates containing: 

a) Average values of damping parameters estimated from three sets of data each with a single frequency, 

C T Nq = 5.9 ± 0.42, C T m =-18.7 ±1.2 

b) Estimated parameters from a combined set of data with all three frequencies included, 

C T Nq =5.5±0.18,C^ =-17.6±0.53 

Because of statistically insignificant differences in the two sets considered it was decided to use the combined data 
including all three frequencies. The final estimates and standard errors for five different test conditions are 
summarized in Table I. 


Table I. Estimated damping parameters from wind tunnel data for BT configuration 

a 0 

(deg) 

a A 

(deg) 

c N 

c 

'-'m 

r T 

L Nq 

<7 

R 2 

c T 

^ mq 

(J 

R 2 

5 

15 

5.5 

0.18 

0.86 

-17.6 

0.53 

0.99 

15 

15 

5.4 

0.22 

0.78 

-15.5 

1.30 

0.96 

5 

10 

5.4 

0.23 

0.70 

-15.4 

0.60 

0.99 

10 

10 

4.7 

0.26 

0.83 

-13.7 

0.85 

0.99 

20 

10 

4.7 

0.37 

0.30 

-16.6 

0.85 

0.86 


Damping parameters in Table I are compared to the out-of-phase components of small amplitude (a A = 3°) 
oscillatory data presented in Figs. 10 and 1 1 . These components are expressed as 

C aq = C aq +C ad a = N or m (23) 

where C N6[ and C m(i are constants for data without frequency effect or they equate to their unsteady counterparts 
when frequency dependence is present 5 . In the given experiment, frequency dependence is noticeable for C aq where 

ct 0 > 10°. 

Assuming that the tail is the main contributor to damping, the out-of-phase components can be formulated 13 as 
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a = N ov m 


(24) 


— 3s 

C a =C a (1 + — ) 

a q V Qa 


From Fig. 4, dslda ~ 0.5, is representative of the slope of s(a) near a 0 = 0. Furthermore, from Fig. 10 and 11, 

C Nq =12.9 and C = -37.9, for -4 ° < a < 4° and all frequencies. Therefore from Eq. (24) the damping parameters 

can be obtained as C N = 8.6 and C mq = -25.3. Both values are higher than those in Table I. The reason for these 

differences could be due to a comparison of results from small and large amplitude oscillatory data and partly due to 
Table I values not containing any contribution to damping from the wing. 



Figure 10. Estimated normal-force damping 
coefficients due to pitch rate from oscillatory 
data for BT and BWT configurations. 



Figure 11. Estimated pitch damping coefficients 
from oscillatory data for BT and BWT 
configurations. 


Rigorous estimation of the out-of-phase components at different ao can be obtained from small amplitude 
oscillatory data at more than five carefully selected frequencies 5 . These results can provide first indication of the 
damping magnitude and its variation with angle of attack. 

B. BW Configuration 

The BW model equations follow from the general model in Eqs. (10-15) as 

C N (t) = C b n (oo; a) + Cx (t;a)+ C™ (oo; a) ^ (25) 


C m (t) = C B n (r Xl -a) + cZ 0 +C%{t,a)x w (a) + C™{K-, a )^ (26) 


As in the complete model, the normal force coefficient was partitioned as 

r w _ r w r w 

L JV - L NA + L NS 

then a time delay due to flow separation can be estimated from 


t n (cc)Cns + Cns ~ Cns (°°; a ) 


(27) 


(28) 
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Similarly, for identification of the pitching-moment model, the coordinate of the wing center of pressure can be 
estimated from 


r x (a)x w +x w = x w (oo;a) (29) 

Now the estimation problem for the C N model parameters can be formulated as follows: given time histories of 
C N (t) , a(t) , and a(t) , and analytical functions C% s ( go ; a) , cj^ A (oo;a ) , and C^(oo;a), estimate unknown 

constants t n and Cy^ . Using simplified notation, r} = C^ s , 0 l =r N , 8 2 = Cy^ , u - Cj^(oo;a) , 

y = C N (t) - (oo; a ) - C% A (oo; a ) , then 

8\>) + r] = u (30) 

y = j]+0 2 d (31) 

A similar formulation can be made for C m where the unknown parameters to be estimated are t x and . The 

unknown parameters in the state-space models were estimated by an output error method explained in Ref. 4. 

Parameter estimates and their standard errors are given in Table II. Both time constant and damping parameters 
were considered independent of a for each data set analyzed. The average value of r v ~ 1 3 for all five cases agrees 
with limited previous experience. The estimates of the second time constant exhibit an increase with increased angle 
of attack. These values could not be verified because no a priori information on them existed. 


Table II. Estimated damping parameters and non-dimensional time constants 
from large-amplitude oscillatory data, BW configuration. 

a n , deg 

a A , deg 

t n 

r BW 
L Nq 

U 

r BW 



0 (a) 

0 (a) 

8 (<T) 

0 (<T) 

5 

15 

12.5 (1.2) 

4.8 (2.2) 

50.0 (3.2) 

5.1 (.63) 

15 

15 

9.6 (1.4) 

8.8 (3.7) 

40.6 (2.6) 

7.4 (.78) 

5 

10 

14.4 (1.9) 

5.8 (2.9) 

25.6 (1.7) 

3.8 (.71) 

10 

10 

12.7 (1.6) 

10.5 (4.0) 

34.0 (1.6) 

10.5 (.69) 

20 

10 

16.6 (1.4) 

13.7 (4.4) 

56.2 (.67) 

0.2 (.67) 


Damping parameters in Table II were compared in Figs. 12 and 13 with the estimates from small amplitude 
oscillatory data. For the low a region, the out-of-phase components are approximately C™ = 5 and C™' = -3 . 

This implies that Cy^ estimates in Table II are very high and that C^f estimates are incorrect due to their non- 
physical values. The reason for these discrepancies was found in the low sensitivity of responses C N (t) and C m (t) 
to both damping parameters. This insensitivity was demonstrated computing these responses and their residuals 
using estimated r N and r x values along with appropriately specified values of Cy^ and . The results, in each 
case, indicated almost no change in the responses and residuals. 
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Figure 12. Estimated normal-force damping 
coefficients due to pitch rate from oscillatory 
data for BW configuration. 



Figure 13. Estimated pitch damping coefficients 
from oscillatory data for BW configuration. 


Measured and computed time histories of C N and C,„ using parameter values from Table II are plotted in Figs. 14 
and 15. Although forced-oscillation experiments were performed separately for each frequency, time history plots 
show all 3 frequency cases stacked into one plot for analysis. In general, there is good agreement between these 
plots. Some effect of uncorrected model error in C„, is visible in part of the data with /= 0.4Hz. Figs. 16 and 17 
show both coefficients plotted against a with /= 0.8 Hz. Included in both figures are also steady values of both 
coefficients. Part of the C v (a) curve for 10° < a < 20° with pronounced unsteady effect is captured quite well by 
identified model despite inaccuracies in both damping parameters. Unexplained differences, however, are apparent 
in the linear part of C v (a) . The agreement between measured and computed C m (a) is poor. Finding an 
explanation for this requires additional study. 



Figure 14. Measured and computed time 
histories of normal force for BW configuration, 
a 0 = 5°, and a A = 15°. 



Figure 15. Measured and computed time 
histories of pitching moment for BW 
configuration, a 0 =5°, and a A = 15°. 
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a 0 , deg 

Figure 16. Variation of measured and computed 
normal force with a for BW configuration, 
a 0 = 5°, and a A = 15°. 



Figure 17. Variation of measured and computed 
pitching moment with a for BW configuration, 
a 0 = 5°, and a A = 15°. 


C. BWT Configuration 

As in the preceding section, model Eqs. (10-11) can be written with simplified notation giving 


d x T] + rj = u 


fl 



_i_ 

^2 

yi_ 



i 

A. 


where 

^1 - T s 0 2 - C N q - C,n q 

rj-s u- s(co) 

h 1 (7j) = S / C N (t;a t ) 
h 2 (7/) = -£ t S t C N (t;a t ) 


(32) 


(33) 


and outputs are defined as 


y ] = C N (t\ a) - C,v (go; a) - Cy (t; a) 

j 2 = C m ( t;a )- Cf (oo; a) - (t; a)x w (t; a) 

Nonlinear functions /?i(? 7 ) and h\ ( // ) are computed by integrating differential equations, Eqs. (28-29), using values 
of t n and r x estimated from BW forced oscillation data. As in the previous cases for this study, it was assumed 
that the unknown parameters are independent of a. 

Unknown parameters, 0 = \r (: C N ■ C m ], were estimated by stacking three sets of oscillatory data for each 

case defined by [ a n , a A ]. Each set in the stack represented only one cycle of data at a single frequency. From the 
stacked data, v e was estimated with an average value of 27 and a relatively high standard error of 5.6. Estimates 
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for C N varied from -5 to 4, and for C n 
values for C, 


values ranged from 14 to 24. Negative values for C N and positive 
l, q 

are unacceptable as follows from the previous discussions and physics of the problem. 

Figures 18 and 19 compare measured and computed coefficients C N (t) and C m {t). The measured data were 
taken from three cycles at a 0 = 20° and a A = 10°. The computed coefficients were based on parameter estimates 
t £ = 26, C Nq = -5, and C mq = 14. The agreement between measured and computed coefficients is very good. 




Figure 18. Measured and computed time 
histories of normal force for BWT configuration, 
a 0 = 20 °, and a A = 10 °. 


Figure 19. Measured and computed time 
histories of pitching moment for BWT 
configuration, a 0 = 20 °, and a A = 10 °. 



10 15 20 25 30 35 10 15 20 25 30 35 



Figure 20. Residuals from measured and 
computed normal-force coefficient for BWT 
configuration, £* 0 = 20 °, and a A = 10 °. 


Figure 21. Residuals from measured and 
computed pitching moment coefficient for BWT 
configuration, £* 0 = 20 °, and a A = 10 °. 


The residuals plotted in Figs. 20 and 21 against a indicate deterministic errors in the data unexplained by the 
model equations. Ideally residuals should be zero mean and normally distributed. Strong insensitivity between 
model damping parameters, and total force and moment coefficient was also found in this case. 

After parameter estimation for BWT configuration, an attempt was made to assess the prediction capabilities of 
the final model. The parameters in this model were selected as shown in the following table. 
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Parameter Specifications for Simulation 

Configuration 

ITT 

II 

C mq =~ 15.2 

BT 

t n -13 

% = 26 

BW 

T e ~ 27 


BWT 


The coefficients C N (t) and C m (t ) were computed for three frequencies, at u {] = 20° and a A = 10°, and compared 
with measured data. As in the previous case where C N (t) and C m (t) were estimated under the same conditions, the 
agreement between prediction and measurement was very good and the residuals were similar to those in Figs. 20 
and 2 1 . 

It is important to note that accurate assessment of model prediction must be based on data sets not used in 
parameter estimation. Validation data sets were not available in this case. 

V. Simulation Results 

Simulation results were used first, to ensure that all estimation algorithms were working properly and second, to 
assess a few experiment design issues that likely produced the strong insensitivity between parameter estimates and 
measurement data in this study. For brevity, only results from using the BW model given in Eqs. (25-29) and a 
limited number of test conditions are provided as examples. Test conditions demonstrate the effect of measurement 
noise, number of oscillation cycles, and model structure error in the form of measurement bias. A number of other 
experiment design factors are important for system identification of nonlinear unsteady systems but these are not 
considered in this paper. Discussion of frequency selection and other test design issues can be found in Refs. 4 and 
14. 

Model parameters were chosen to reflect a realistic set of values that are supported by the model identified in this 
paper. All simulation results are shown for a forced-oscillation experiment at do = 15 degrees and aA — 15 degrees. 
Two polynomial cases are demonstrated for the non-dimensional time constant, r N , shown in Table III. In all 
simulation examples the starting values for parameter estimation are chosen to be 10% in error. This is a small 
starting error, but it is sufficient to demonstrate the estimation algorithms and specific test conditions of interest. 


Table III. Simulation model parameters applied to BW configuration at a 0 = 15 deg and a A = 15 deg. 

Model Parameters 

Oscillation 

Frequency 

r N 

r BW 
L Nq 

Tt 

r BW 

mq 

Parameter Structure 

/(Hz) 

c 0 + c x a 

c 0 

c 0 

c 0 

Case 1 (constant) 

0.8 

15 + 0a 

0.5 

20 

-2 

Case 2 (linear poly) 

1.6 

5 + 19. 0986a 

0.5 

20 

-2 


Case 1, with all parameters equal to constants, is an example demonstrating the efficacy of the estimation algorithms 
and the relative magnitude of the steady flow damping term, C ; (/ . Figure 22 shows the results of using simulated 

measured data to estimate the normal force coefficient and predict its response against angle of attack. Model 
parameters were estimated exactly using only 1 oscillation cycle of data at 0.8 Hz, providing some confidence in the 
estimation algorithms. 

Because the steady flow damping coefficient, C ™ , has a small value the damping term contributes relatively 

little to the response and consequently requires careful testing for successful model identification. For comparison. 
Fig. 23 shows the same model prediction as in Fig. 22 except only the steady flow damping term is present. The 
difference of the predicted response from the static curve is not easily visible demonstrating the very small 

contribution of the steady flow damping term. If is the sole term used to explain the response during model 

identification a very large and unrealistic value results indicating that this particular response can only be explained 
by a model with unsteady dynamics. 
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The predicted response shown in Fig. 23 uses measured data for initial conditions to solve BW Eqns. (25-29), 
consequently the predicted response exhibits a small transient at a = 15 degrees. This response reflects the transition 
from the measured initial condition to the predicted steady harmonic motion. 



Figure 22. Variation of measured and computed Figure 23. Variation of measured and computed 

normal force for BW configuration, at a 0 =15°, normal force without unsteady model component 

and a A = 15°. for BW configuration, at a 0 = 15°, and a A = 15°. 


To increase the realism of the estimation problem, possibly reflecting a change from steady to unsteady regions 
between low and high angle of attack, the unsteady time constant is allowed to vary with angle of attack. This is 
labeled as Case 2 in Table III. In this case the time constant is defined linearly from 5 to 15 over a range of angle of 
attack from 0° to 30°. Initially only 1 cycle of oscillation data is used to simulate the real data case investigated in 
this paper. Two runs are made, one without noise and one with a very modest level of noise so that the signal to 
noise ratio is approximately 60. Measured and computed normal force coefficients are shown in Figs. 24 and 25 for 
this example. In both cases the estimated model predicts the response well. In the case with noise the parameters are 
estimated with extremely poor accuracy as shown in Table IV. 


Table IV. Estimated damping parameter and non-dimensional time constants 
from simulated large-amplitude, 1-cycle oscillatory data, BW configuration at/=1.6Hz, a Q = 15, a A = 15 deg. 

Parameter true values: = 0.5, c 0 = 5, q = 19.0986, where r A r=c 0 +c 1 a 

Signal/noise 

F it error 

# cycles 

Aj 

II 

<9 = c 0 

9 = q 




e (a) 

e {&) 

9 (a) 

OO 

3.9e-06 

1 

0.4997 (0.0004) 

4.9995 (0.0004) 

19.1001 (0.0010) 

60 

0.017466 

1 

1.8424 (1.8094) 

4.9343 (1.6240) 

15.1883 (4.6112) 


The modest amount of noise masked the limited information content of the data and caused large parameter 
insensitivity. This is especially a problem for the very small steady flow damping term, , in this case. To 

improve the model estimation in the presence of noise, the number of cycles of oscillation analyzed is increased and 
varied from 6 to 14 cycles. The parameter estimation results, shown in Table V, demonstrate, as expected, that 
increasing the number of oscillation cycles improves both parameter estimates and their standard errors. However, 
unexpectedly, the number of cycles required to improve parameter accuracy is relatively large compared to 
conventional practices. 
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Off deg 

Figure 24. Variation of measured and computed 
normal force for BW configuration, at ar fl =15°, 
and a A = 15°. Normal force time constant varies 
linearly with angle of attack. No noise is present. 



Figure 25. Variation of measured and computed 
normal force for BW configuration, at a,r 1 5°, 
and a A = 15°. Normal force time constant varies 
linearly with angle of attack. Signal/noise ratio 
=60. 


Table V. Estimated damping parameter and non-dimensional time constants 
from simulated large-amplitude oscillatory data, BW configuration at f=l. 6Hz, « 0 =15, a A = 15 deg. 

Parameter true values: C™ = 0.5, c 0 = 5, c, = 19.0986, where r N = c 0 +q a 

Signal/noise 

F it error 

# cycles 

^ 13 , 

II 

II 

o 

6 = cj 




0 (a) 

6 (a) 

e (a) 

00 

1.0e-06 

6 

0.50 (0.0001) 

4.999 (0.0001) 

19.099 (0.0002) 

60 

0.0158 

6 

0.3 (0.68) 

4.5 (0.62) 

20.7 (1.75) 

60 

0.0158 

9 

0.6 (0.55) 

4.6 (0.50) 

20.0 (1.43) 

60 

0.0155 

14 

0.5 (0.43) 

4.8 (0.40) 

19.7 (1.12) 


It is common practice in the aerodynamics field to apply one-factor at a time (OF AT) procedures when 
performing wind tunnel experiments. For example, typically, angle of attack is adjusted systematically over some 
range of interest. This lack of randomization can introduce systematic or bias errors into an experiment. To assess 
this issue simulated data was generated for 2 values of measurement bias error on the normal force coefficient. The 
bias errors were computed as a percentage of the maximum value of the normal force coefficient time history and 
the percentages were selected as 0.1% and 1%. Table VI shows the estimation results when bias error is included for 
the same simulation considered in Table V. For each case with bias, 14 cycles of oscillation are used in the analysis. 
Estimation results using a model with bias error but without measurement noise are shown in the first three rows of 
Table VI. The results are revealing and indicate that even a very modest bias error can degrade estimation accuracy 
while a bias of 1% of maximum signal value confounds the estimation process, producing a non-physical mean 
value for the steady flow damping term. Estimates in the last 3 rows show that the modest bias error in combination 
with measurement noise significantly degrades the estimation accuracy and the 1% bias completely defeats the 
estimation process by preventing any convergence of the estimation algorithm. These results are similar to problems 
encountered with the real wind tunnel data in this study. 
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Table VI. Estimated damping parameter and non-dimensional time constants 
from simulated large-amplitude oscillatory data, BW configuration aty=1.6Hz, «o=15, a A = 15 deg. 
Parameter true values: C^' = 0.5, c 0 = 5, q = 19.0986, where t n = c 0 + c , a 

Signal/noise 

Bias error 

Fit error 

II 

0 -c 0 

9 = q 


(% max) 


9 (<r) 

6 (a) 

9 (a) 

00 

0 

1.0e-06 

0.50 (0.0001) 

4.999 (0.0001) 

19.099 (0.0002) 

00 

0.1 

0.0013 

0.46 (0.036) 

4.8 (0.03) 

19.8 (0.094) 

00 

1.0 

0.0129 

-0.1 (0.36) 

2.9 (0.33) 

26.5 (0.94) 

60 

0 

0.0155 

0.5 (0.43) 

4.8 (0.40) 

19.7 (1.12) 

60 

0.1 

0.0154 

0.4 (0.43) 

4.7 (0.40) 

20.6 (1.11) 

60 

1.0 

~ 

~ 

~ 

~ 

indicates estimation algorithm did not converge. 


VI. Concluding Remarks 

This paper first examines aircraft model identification from wind tunnel data of an airliner and its components: 
body, body-tail, body-wing, and the complete model. A mathematical model of an aircraft with structure sufficiently 
general to model nonlinear unsteady behaviors is considered. The measurement data included results from 
conventional static and one degree-of-freedom in pitch oscillatory tests over a wide range of angle of attack but 
limited frequencies and amplitudes. 

All models with estimated parameters fit the measured data quite well. Damping parameters estimated from BT 
large-amplitude oscillation data were consistent with small-amplitude out-of-phase measurements and after 
including an a-correlated term to correct for an inconsistency between normal force and pitch moment 
measurements, a very large percentage of the variation in the data was explained by the model. Damping parameters 
for the BW and BWT configurations were poor, however, the values were explained by a lack of parameter 
sensitivity or information content of the data. 

Conventional forced-oscillation experiments generally focus on obtaining averaged single-cycle data at a limited 
number of frequencies. This approach facilitates a straightforward (sometimes graphical) approach for estimating 
out-of-phase damping in steady-flow linear regions. Identification of more complex models that capture nonlinear 
unsteady behaviors is not well suited for use with conventionally designed experiments as demonstrated by the very 
limited success of the BW and BWT model identification in this study. Simulated data were used to confirm the lack 
of parameter sensitivity in a conventional experiment and to demonstrate the potential severity of estimation error 
due to measurement noise and bias error. Although, theoretically, certain harmonic information is fully contained in 
a single cycle of data, both the actual and simulated results of this study support the need for a much greater number 
of cycles of data to overcome the effect of measurement error. This is especially true for situations where it is 
necessary to estimate model terms with relatively small magnitudes compared to the total aerodynamic force or 
moment. 
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